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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 
RESEARCH MEMORANDUM 

CALCULATIONS OF LAMINAR HEAT TRANSFER AROUND CYLINDERS OF ARBITRARY 

! 

CROSS SECTION AND TRANSPIRATION-COOLED WALLS WITH 
APPLICATION TO TURBINE BLADE COOLING 
By E. R. G. Eckert and. John N. B. Livingood 


SUMMARY 

An approximate method for the development of flow and thermal 
"boundary layers in the laminar regime on cylinders with arbitrary cross 
section and transpiration- cooled walls is obtained by the use of von 
Harman’ s integrated momentum equation and an analogous heat -flow equa- 
tion. Incompressible flow with constant property values throughout the 
boundary layer is assumed. The velocity and temperature profiles within 
the boundary layer are approximated by expressions composed of trigo- 
nometric functions. Shape parameters for these profiles and functions 
necessary for the solution of the boundary-layer equations are presented 
as graphs so that the calculation for any specific case is reduced to 
the solution of two first-order differential equations. 

The method is applied to determine local heat-transfer coefficients 
and surface temperatures in the laminar flow region of transpiration- 
cooled turbine blades for a given coolant flow rate, or 'to calculate the 
coolant flow distribution which is necessary in order to keep the blade 
temperature •uniform along the surface. 


INTRODUCTION 

Transpiration cooling is a very effective means for keeping sur- 
faces which axe subject to a hot gas stream at a low temperature. For 
use of this method, the surface is fabricated from a porous material 
and a cooling fluid is blown through the pores. Along the outside sur- 
face the cooling fluid builds a film which insulates the wall from the 
hot gas stream. The transpiration-cooling method may be applied to the 
cooling of structural parts in propulsion systems such as gas-turbine 
blades, combustion-chamber walls, and rocket nozzles. If a heated fluid 
is blown through the porous wall, -the same method may be used to keep 
the surface tenperature of the wall at a value which is higher than the 
temperature in the outside flow. In this way, it may be applied for 
de-icing of wings and other parts of airplanes and in the propulsion 
system. 
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This report presents a method by which the heat transfer connected ' 
with transpiration cooling in two-dimensional 1 ami nar flow around bodies 
of arbitrary cross section can be calculated. It considers only the 
case where the fluid blown through the porous wall is the same as the 
one in the outside flow. Although the particular application considered 
in this report is the cooling of turbine -blades, the method itself may 
be applied to other applications as well. The procedure by which the 
calculation may be carried out for any particular application is 
described in an appendix in this report. 

The determination of the heat transfer is based on the calculation 
of the thermal boundary layer which builds up around any body in a flow 
field. This thermal boundary layer is interconnected with the flow 
boundary layer for variable fluid properties or superimposed on the 
velocity boundary layer when the properties are independent of tempera- 
ture. The build-up of the flow boundary layer is determined by the 
pressure distribution around the body under consideration. For the type 
of pressure distribution occurring in the aforementioned applications, 
only approximate methods of calculation are straightforward enough for 
engineering purposes. For the determination of the thermal boundary 
layer in particular, two types of approach are known. 

The first approach was introduced by Krouj iline for the calculation * 

of heat transfer on solid surfaces and presented in more detail by Fros- 
ling and others (reference l) . In this method, the flow boundary layer 
has to be known before the thermal boundary layer can be calculated. 

Usually the method, as introduced by von Karman (reference 2) and 
Pohlhausen (reference 3) , which fulfills the integrated momentum equa- 
tion of the boundary layer, is applied for the calculation of the flow 
boundary layer. More recently, Wieghardt and Walz (reference 4) have 
used, in addition to the momentum equation, an integrated energy equa- 
tion and Tetervin and Lin (reference 5) have introduced a still more 
general integral condition for the boundary layer which may be used in 
such calculations. The use of these expressions gives better agreement 
with exact calculations and with measurements in special cases, par- 
ticularly in regions where the pressure increases in flow direction. 

Since, in the application' considered, the regions- of most concern are 
those where the pressure decreases, the integrated momentum equation, 
which Is simplest to handle, will be used. Schlichting (reference 6) 
used this equation to calculate the flow boundary layer on a porous 
surface through which fluid is sucked in order to keep the boundary 
layer l aminar or to prevent flow separation. A paper by Dorodnitzyn 
(reference 7) extended the method to include the effect of Mach number 
and of variable property values. However, this extension is developed 
only for zero heat transfer. The calculation of the flow boundary 
layer in this report will essentially follow Schlichting* s approach. 


2226 


NACA EM E51P22 


3 


After the flow boundary layer is determined, the thermal boundary 
layer can be calculated according to the method of Kroujiline by use of 
an integrated heat-flow equation. Velocity profiles known from the cal- 
culation of the flow boundary layer and temperature profiles within the 
boundary layer whose shapes are approximated by a polynomial expression 
are introduced into this heat-flow equation and the equation is solved 
for the thickness ratio of the thermal boundary layer to the flow bound- 
ary layer (reference l) . This procedure, however, becomes quite 
lengthy, especi ally when the number of terms in the polynomial expres- 
sion for the temperature profile .is increased in order to improve the 
accuracy of the approximation. When an attempt was made at the HACA 
Lewis laboratory to extend this method, which was previously used only 
to deter min e the heat transfer on solid surfaces, to transpiration- 
cooled porous walls, it was found that the procedure for the solution 
of the heat-flow equation becomes much simpler and more straightforward 
if the equation is solved for a thermal boundary-layer thickness called 
convection thickness. This will be explained in detail in a later 
section of this report. After the method was developed and the calcula- 
tions were finished, a note by Dienemann (reference 8) was found in 
which the same kind of approach is briefly described. Dienemann applies 
the method to calculate heat-transfer coefficients on solid surfaces 
and proposes to extend it in such a way that a temperature variation 
along the solid surface and the influence of internal frictional heating 
within the boundary layer may be accounted for. He also shows that this 
method is superior to other approximations, which were compared by 
Goland (reference 9) with an exact solution for a cylinder with circular 
cross section. It may be of interest to mention that Goland obtained 
the exact solution from the fact that the differential equation describ- 
ing the temperature boundary layer around an infinite cylinder in a flow 
normal to its axis and for a fluid with a Prandtl number of 1 has 
exactly the same form as the differential equation describing the span- 
wise flow within the boundary layer on a yawed infinite cylinder. Cal- 
culations of the heat transfer on a transpiration-cooled flat plate 
which considered the variation of property values with temperature were 
made by Yuan (reference 10) with the assumption that the total- 
temperature profile within the boundary layer is similar to the velocity 
profile. This assumption is valid for a Prandtl number equal to unity 
and for no pressure gradient. The purpose of the present investigation 
is to consider the influence. of pressure gradients as well as Prandtl 
numbers different from unity. 

A second approach for obtaining an approximate solution of the 
thermal boundary layer was described in reference 1. It is still 
simpler than the approach by Kroujiline, since in this method it is not 
necessary to calculate the flow boundary layer prior to the determina- 
tion of the thermal boundary layer. This method uses exact solutions 
of the boundary -layer equations which are known- for a special type of - 
pressure variation in the flow which is encountered on wedge-shaped - •- 
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"bodies. These velocity profiles and the temperature profiles are used 
to approximate the actual profiles for arbitrary pressure variations. 

A differential equation is set up with the condition that the growth of 
the boundary layer on any place of the cylinder with arbitrary cross 
section investigated be the same as for the wedge-type flow when the 
boundary-layer thickness and the pressure gradient have the same values 
in both cases. When this idea is applied to the momentum thickness of 
the boundary layer, the resulting equation is identical with Harman' s 
integrated momentum equation. When it is applied to the convection 
thickness of the thermal boundary layer, the resulting differential 
equation fulfills the condition that the heat transferred from the sur- 
face to the fluid has to be found again within the boundary layer 
(reference ll) . This method was compared with exact solutions and 
experimental values in reference 1 and in investigations performed at 
the University of California (reference 12), and the agreement obtained 
was quite satisfactory. This method can easily be extended in such a 
way that it includes effects of variations of the surface temperature 
and of internal heating (reference ll). For such an extension, the 
corresponding exact solutions for the wedge-type flow have to be known. 
Such solutions, which take into account the effect of a surface tem- 
perature variation and of internal heating, are presented in refer- 
ences 11 and 13. The method may also be extended to the transpiration 
cooling of porous surfaces, as soon as the corresponding exact solutions 
for this case are known. A few such solutions are presented in refer- 
ence 14; However, too few solutions are given for use as a basis for 
the approximate method. Brown (reference 15) has recently made an 
extensive calculation to obtain exact solutions in transpiration-cooled 
porous surfaces of the wedge-flow type which include the effect of 
pressure gradients and of variable property values. The results of this 
calculation are now being used to extend the method mentioned in the 
preceding paragraph to transpiration cooling. 

This paper deals with the method described as the first type of 
approach. The method has the advantage that it may be applied to cases 
for which the corresponding wedge-type flaw and heat transfer are not 
known. 

STATEMENT OF PROBLEM AND SIMPLIFYING ASSUMPTIONS 

This report is a contribution to the problem of determining the 
development of the thermal boundary layer and the heat-transfer coeffi- 
cients on a body of arbitrary cross section with porous walls in a two- 
dimensional flow. Figure 1 shows the cross section of such a body. At 
Reynolds numbers which are sufficiently high, the flow around such a 
body may be subdivided into the boundary-layer region which surrounds 
the body with a very small thickness, and the potential flow, which 
determines' the pressure distribution around the body. The highest 
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pressure on the body is found at the stagnation point. The pressure 
then decreases in flow direction along both sides of the body, and 
usually increases again, later on. The pressure variation along the 
body deter min es the development of the flow boundary layer and also 
whether the boundary layer is laminar or turbulent. Usually, the lami- 
nar part is confined to a region near the nose of the body. The region 
with laminar boundary layer is investigated in this report. The flow of 
coolant through the porous surface may be described by the velocity v w 
with which the coolant leaves the surface in a normal direction. The 
questions which will be answered in this report are: First, what are 

the local heat-transfer coefficient and the surface temperature on any 
point along the body for any prescribed distribution of the coolant 
velocity v w ; and second, what distribution of the coolant velocity v w 
gives a desired distribution of the heat-transfer coefficient and of the 
surface temperature around the body. Usually, for instance, a constant 
wall temperature is most desirable and the problem is to find that dis- 
tribution of the coolant velocity v w which results in a constant wall 
temperature . 

A number of simplifications have to be introduced in order to keep 
the time required for the solution of a special problem down to a 
tolerable amount. The assumptions made are: The flow is two- 

dimensional and in steady state, internal, frictional heating within the 
boundary layer can be neglected, and property values (density, vis- 
cosity, and heat conductivity) may be considered constant. The influ- 
ence of temperature gradients along the surface of the body is neg- 
lected. Although this influence may be considerable (reference 13), 
there is no quantitative information available for transpiration-cooled 
surfaces which would permit this effect to be taken into account. 

In applying the method to the determination of the temperatures of 
transpiration-cooled turbine blades, neglecting internal frictional 
heating should be admissible, since the temperature differences within 
the boundary layer generated by the cooling process are considerably 
larger than those generated by aerodynamic heating. This will be shown 
in more detail later on. On the other hand, the large temperature 
variation m the boundary layer is connected with a considerable varia- 
tion of the property values. The influence of this variation on the 
heat transfer may be approximately corrected by use of the results in 
reference 15. 


BOUNDARY -LAYER EQUATIONS 

In a coordinate system, shown in figure 1, where the x-axis runs 
along the surface and the y-direction is normal to the surface and under 
the assumptions mentioned in the previous section, the differential 
equations describing the velocity and the temperature within the 
boundary layer are 
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pu* 



pv* 


Su* 

3y* 



Sfu* 

ay * 2 


( 1 ) 


PCpU* 

(All symbols are defined m appendix A.) The equations have to ful- 
fill the following boundary conditions. At the surface (y = 0 ) the 
velocity component u* parallel to the surface is 0, the velocity 
component v* normal to the surface has a finite value v w *, and the 
temperature of the surface is interconnected with the velocity v w * by 
the over - all heat-transfer process so that one of the two values may be 
1 prescribed. At the outer edge of the boundary layer (y =oo) the u* 
component of the velocity transforms asymptotically into the stream 
velocity U* and the temperature transforms into the stream tempera- 
ture. Since only temperature differences appear in equation (3) the 
temperature level does not enter into the problem. Therefore, all tem- 
peratures will be measured from the temperature in the stream as refer- 
ence temperature, and t will be interpreted as the temperature differ- 
ence from this reference temperature. Consequently, outside of the 
boundary layer t = 0. 

' In a gas stream, it is advantageous to interpret t as the total 
temperature difference. In this case, equation (3) already includes the 
effect of the frictional heating for a gas with a Prandtl number of 1. 
Since for all gases, the Prandtl • number does not deviate much from the 
value 1, equation (3) also gives a" good approximation to the real con- 
ditions for gases as long as the temperature differences impressed upon 
the boundary layer by a cooling process are larger than the temperature' 
differences created by internal friction. 

In order to reduce the number of parameters, the differential 
equations will be made dimensionless. For this purpose, all lengths 
measured parallel to x are divided by L, the distance between the 
stagnation point and the trailing edge of the body measured along the 
surface, and all velocities in this direction are divided by an upstream 
velocity Uq. All lengths and all velocities parallel to y are in 
addition multiplied by the square root of Reynolds numbers Reg based 
on the body length L and the upstream velocity Uq: 



Su* Sv* _ n 
Sx* Sy* 


St * St , 

_ + pcpv _ = k 


S 2 t 


dy 


*2 


( 2 ) 

(3) 
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The pressure gradient bp* /bx* is impressed upon the houndary layer by 
the potential flow outside the boundary layer and can be replaced by , 
the stream velocity gradient by use of Bernoulli^ equation 

< 5 > 

In this way, equations (l) to (3) transform into 

p 

5u , du TI 5U , b ^ u 

u S + T s7 = D 5J + i7 (6) 


— + — = 0 
dx by 


bt , bt 1 d 2 t 

u r— + v — = — 

ox dy Pr ^ y 2 

The hound ary conditions for these equations become 

y = 0 u = 0 v = v w (x) t = t^Cx) 


y — > oo u — * U t — * 0 (l( 

The equations are now integrated over y from y = O to y = oo. The 
result of this integration is Kerman's integrated momentum equation 

■&£ (u - u)u ^ + (u - u) ^ - T »° - (^) w (1] 
and the heat-flow equation 


1 /bt 


s/ tudy- 


In order to simplify these equations, the following characteristic 
boundary-layer thicknesses are introduced: 

(l) the displacement thickness of the flow boundary layer 


5 d J (1 - H) ay 
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(2) the momentum or impulse thickness of the flow boundary layer 



(14) 


and (3) the convection thickness of the thermal boundary layer 



(15) 


where 0 is the ratio t/l^ with the limiting values 0 = 1 for 
y = 0 and 0 = 0 for y =». By introduction of these boundary-layer 


thicknesses and with ^ written as U’ , equations (ll) and (12) are 


transformed into 


_d_ 

dx 


(UV + 5 d UU' - v„U 

d / . ttc \ 4- ^w ( bd\ 

ax (V^tjc) - t w v w - - Pr^ y J w 


(16) 

(17) 


VELOCITY AND TEMPERATURE PROFILES 


To obtain an approximate solution of the last two equations, 
approximate expressions for the velocity and temperature profiles will 
be introduced. The accuracy of the results of this calculation will 
depend on how well the actual profiles are approximated by the assumed 
shapes. The temperature profiles as well as the velocity profiles are 
chosen as a one -parametric family. The parameter for each f ami ly is 
determined in such a way that the assumed profiles fulfill the exact 
boundary -layer equations (6) to (8) at the wall surface. 


fbu) _ TT dU ( c^u \ 

Vw W/v = U Sx U y 2/ W 


(18) 


be 




_i_ 

Pr 


a 2 0\ 


dy Z /w 


(19) 


r 

c 

c 

c. 
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In addition, the following boundary conditions will be fulfilled 


II 

0 

0 

11 

2 

cr> 

II 

H 

(20) 

II 

8 

S3 

II 

2 

0 

II 

CD 

(21) 


Originally, polynomial expressions were used to approximate both the 
velocity and the temperature profiles. Schlichting, however, pointed 
out that better approximations may be obtained by expressions composed 
of trigonometric functions (reference 6) . The following profiles are 
used in this report: 

0 < |<1 
- e' 3 ^ - sin f g) 

r ( 22 ) 


j 

By use of this profile, the following expression is found from equa- 
tion (18) for the shape parameter A: 


Velocity profile . - For A > 0, 
and for ^ — 1 

H = 1 - Ae 
U 


U’B 2 - | v w & 
9 + (3-f)v, ; 8 


(23) 


The velocity profile equation (22) has already been used by Schlichting 
(reference 6). For a solid flat plate (u* = 0 and v w = 0) it 
approximates the Blasius velocity profile very well and for U 1 = 0 
and v w = - 3/5 transforms into the exact asymptotic suction profile 
as calculated by Schlichting (reference 6) . For A < 0, 0<^<1 

H = sin — ^ + A sin — 2C ( 1 - sin — 

U 2 & 2 8 \ 2 8/ 


and for g — 1 


> 


u 

U 


1 


J 


(24) 
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The shape parameter in this case. 


as obtained from equation (18), 


is 


A = 


U«5 2 



it * 
2 


* r 
— v^o 

2 w 


(25) 


For the solid flat plate, this profile gives the same expression as 
equation (22). It will be shown later that it approximates the separa- 
tion profile as calculated by Hartree (reference 16) better than the 
usual polynomial expression. Separation occurs at A = - 1 and the 
corresponding profile has the shape 


u . c it y 

U = 6 “ 2 2 I 


(26) 


Temperature profile . - For 0 < < i 


1 - = e 


% 


and for r > 1 


= 1 - 




sin 


£ JL 

2 6 + 




&t 


r (27) 


0=0 


The shape factor for this profile is obtained from equation (19) as 
follows: 


= 

^ it + Pr v w 5 t 

In this case, a negative shape factor A(- is always connected with a 
positive value of v w whereas for the velocity profile, the sign of the 
shape factor depends, in addition, on the pressure gradient dU/dx. 

With these profiles, the different boundary-layer thicknesses and 
the velocity gradient at the wall may be calculated. 
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^ = 0.1366 + 0.03791 A - 0.00786 A 2 J> A > 0 


du\ = U 
dyjw & 


3-fiA 


J ^ 






= 0.1366 - 0.01456 A - 0.02618 A 2 A < 0 


f ^\ -«2(l + A) 

y dy/w 2 5 


J 


The nondimensional temperature gradient at the wall is 


d 0 ] * 1 / n . \ 

dy/w " 2 S t ^ + ^ 


11 

(29) 

(30) 

(31) 

(32) 

(33) 

(34) 

(35) 


The corresponding expressions for the convection thickness of the 
thermal boundary layer are presented in appendix B. 


TRANSFOEMATION OF BOUNDARY -LAYER EQUATIONS 


Flow Boundary Layer 


Multiplication of equation (16) by 
tion of the first term give 


6i 

U 


and a partial differentia- 




dBj 

dx 


2 + 


B d 

6 , 




V v^i “ u 



(36) 


The expression in the parentheses of the second term on the left-hand 
side and the term on the right-hand side are functions of the shape 
parameter A. Therefore, the momentum equation for the boundary layer 
may be written in its final form: 
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I = f 2 ( fi > - f l < A > “V + V*i 


(37) 


with the following expressions for the two functions f-^ and fg which 
are obtained from equations (29) to (34): 


„ /aN „ 0.3634 - 0.03005 A 

fl (A) = 2 + p 

x 0.1366 + 0.03791 A - 0.00786 Ar 




(38) 


>A> 0 


f 2 (A) = [| + ( 3 - §)A] (0.1366 + 0.03791 A - 0.00786 A 2 } 


(39) 


fl (A) = 2 + 


0.3634 - 0.1366 A 




0.1366 - 0.01456 A - 0.02618 A 


f P (A) = J (l + A) (0.1366 - 0.01456 A - 0.02618 A 2 ) 


J 


(40) 

< 0 

' (41) 


These functions are presented in figure 2. All curves have a break at 
A = 0 because different expressions approximate the velocity profiles 
for positive and negative A values. Equation (37) is a, linear first- 
order differential equation from which the momentum thickness of the 
boundary layer can be obtained by integration, as soon as the gradient 
U* of the stream velocity U and the porous flow characterized by v w 
are known as functions of x. In order to make the calculations more 
convenient, the shape parameter A can be expressed as a function of 

the two quantities U'&i and VyS 


A>' 0, A = 


p Sj. 

U*5i 2 - - v^ — 


6 7VT 


9 It 


, n ] - &i 

3 " 2 V * 8 i T 


(42) 


A< 0, A = 


IJ'&i - % Vw 5i -± 


* 2 ( 

V 

) 2 + I vS 

Si 

2 ' 

V 6/ 

1 + 2 v w°i 

5 


(43) 


l / 

The ratio -g- occurring in these equations is a function of A (see 

equations (30) and (33) and fig. 2). The functional relation for A is 
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plotted in figure 3. By use of figures 3 and 2, the integration pro- 
cedure for the differential equation (37) becomes very simple. The 
step-by-step procedure for such a calculation is presented in 
appendix C . 

In order to start the calculation at the stagnation point, it is 
necessary to know the boundary -layer thickness 5^ at this location. 
At the stagnation point, the stream velocity U is 0. Since on a 
blunt nose the increase of the boundary-layer thickness is never infi- 
nite, the term on the right-hand side of equation (37) has to be 0. 
This gives an equation 


f 2 (A) - f ± (A) U'Bi 2 + v w 5 ± = 0 


(44) 


Qn the other hand, figure 3 also applies to the conditions at the stag- 
nation point. From both relations U'5^2, a, f-^, and f 2 were each 

obtained as a function of v-y^ by a trial -and-error process. These 
values are presented in table I. The connection between U’ 8 ^ 2 , A, 
and v w 5^ is also shown in figure 3 as the stagnation line. Table I 
cannot be U 6 ed immediately to start the calculation since the value 

( v w s i ) 2 v w 2 

v^i is not known. However, the value v t = -=ry which conta ins 


U'8-i 2 


only known values, is also a function of U' 8^ 2 and is plotted in fig- 

v 2 

ure 4. From this figure U' 8 ^ can be determined for a given -^7 
and 5^ can then be obtained from this value. 


Thermal Boundary Layer 

With the help of equations (28) and (35), the heat-flow equa- 
tion (17) is transformed into 

ii (D5i ) ! , « 1 (, 

v w dx - i + 2Pr Kvyd y- ~ n + ?r K Vw S j “ v w t w dx 


(45) 


' h 

where K denotes the ratio — . The last term on the right-hand side 

o 

accounts for a variation of the temperature along the surface. A dif- 
ficulty arises in this connection. It is known from the results of 

references 11 and 13 that a temperature gradient also has a 

dx 


I 
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pronounced effect on the shape of the temperature profile. Since no 
such effect was included in the assumed profile (equation (27)), it is 
doubtful what significance the last term has. It was therefore neg- 
lected herein, restricting this report to cases where the variation of 
the surface temperature is kept small either by internal conduction 
within the well or by proper choice of v w . The investigation of the 
influence of large surface temperature gradients will be left to future 
work. With this simplification equation (45) can be written in its 
final form 


S (TOt.c) = V 3 (Pr Kv w &) 


(46) 


where f 3 is a function of the product Pr Kv w 8, which is presented in 
figure 5. For a solid wall (v w = 0) the equation simplifies to 

dx = 2 K Pr 5 ^ 


The shape parameter A^ may be written in the form 


At = 


- Pr KvyS 
jt + Pr Kv w 6 


(48) 


This equation, together with the expressions for 


6 . 


hui 

8 


(appendix B) , 


U. 

determines a functional relation between K, ^ , A, and Pr which 


is presented in figure 6. The figure presents the ratios 


J t,c 

8 


for 


three values of the shape parameter A (-1, 0, and +l) . For inter- 
mediate values of A linear interpolation in the range A = -1 and 
A = 0, or A = 0 and A = 1, at a constant value of K may be used 
with good accuracy. The heat-flow equation (46) is again a linear 
first-order differential equation from which the thermal convection 
thickness is obtained by integration when the stream velocity U, the 
flow through the porous surface v w , and the Prandtl number Pr are 
prescribed, and the boundary layer thickness 8 and the shape factor 
of the flow boundary layer A are known from a preceding solution of 
equation (37). The step-by-step procedure for such a calculation is 
explained in appendix C. 


In order to start the calculation, the convection thickness &t,c 

must be known at the stagnation point. A partial differentiation of the 
left-hand side of equation (46) gives 
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d6 t c 

u -£r + u,& t,c - V3 


At the stagnation, point, U is 0, and the boundary-layer increase 

t>c is not infinite on a blunt-nosed body. This results in the 
dx 

equation 


or 


EL 



f 3 (Pr Kv w 6) 


U»8 2 frt.c 
v v 5 6 


f 3 (Pr Kv w 6) 


(49) 


This is a relation between the two unknown ratios 
second relation is given by figure 6. 


and 


't,c 
6 

From both, the values 


K. A 

hzl 

8 


and K can be determined by a trial-and-error procedure, 
are presented in table II and figure 7. 


The values 


When the thermal convection boundary-layer thickness is known, the 
local heat-transfer coefficient follows by a simple calculation. The 
equation which defines the local heat-transfer coefficient h is 


htw “ k, §?) v 


(50) 


This gives for the local Nusselt number, based on the length L, 


..^55 fc£\ 


hL = _ 

^ "W ^ dy*/ 


w 


*/ w 


(51) 


The introduction of equation (17) for the temperature gradient leads to 
an expression 

= ** “ v w] = v w ** ( f 3 “ x > 


(52) 
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Without sweat cooling, the equation simplifies to 


(53) 


In many cases, the whole amount of heat transferred from the out- 
side flow to the wall surface is picked up by the cooling fluid on its 
flow through the porous wall. In this case, the wall surface tempera- 
ture may be calculated by the equation 


Nu it 1 

/R ^ =2K5 


htw — pCpVv. (t c ~ ty-) 


(54) 


where t c is the temperature with which the coolant enters the porous 
wall. This gives for $, the ratio of the difference in gas tempera- 
ture minus wall surface temperature to the difference in gas tempera- 
ture minu s coolant temperature, the expression 




(55) 


COMPARISON OF ASSUMED VELOCITY AND TEMPERATURE 
PROFILES WITH EXACT SOLUTIONS 

The accuracy of the method depends on how well the assumed profiles 
approximate the actual ones. It is therefore necessary to check the 
degree of this approximation with results of exact calculations to the 
extent that these are available. Such a comparison will be made in this ■> 

section. 


' Velocity Profiles 

For the solid surface (v w =0), the comparison can be made with 
exact solutions which were obtained by Hartree (reference 16) for wedge- 
type flow for which the stream velocity is a power function of the dis- 
tance from the stagnation point (U = U^x®) . A set of such profiles 
taken from reference 16 is presented as a family of dashed lines in 
figure 8 where the ratio of the velocity u in the boundary layer to 
the stream velocity U outside of the boundary layer is plotted over 
the dimensionless distance 
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z 




and the shape parameter p is defined "by the equation 
expression for the dimensionless displacement thickness is 

This gives 


a 2m 

P = 5?T 


The 


- &1 [H 

5l 


Pz 


2 

l 



2 U 
x 


A differentiation of U = U^x 111 gives 


TTI _ m U _ 3 U 


Therefore, 


p Z;L 2 = U'5 ± 2 (56) 

The value z ^ as a function of p is presented in reference 1. Since 

o 

U'5^ is a function of A, according to equations (42) and (43), equa- 
tion (56) presents a relation between p and A from which A may be 
calculated for any value of p. The shape of the approximate velocity 
profile is determined by A (equations (22) and (24)). In this way, 
the solid profiles of figure 8 were determined. They axe superimposed 
on the exact profiles in such a way that the displacement thickness is 
the same in both cases. Figure 8 shows that generally the agreement 
between the exact profiles and the approximation used in this report is 
satisfactory. Only for the separation profile (p = -0.1988) are the 
deviations larger, but even for the separation profile the approxima- 
tion by equation (24) is better than the usual four-term polynomial 
approximation, which is shown as a dash-dot line. 

For a porous wall, some exact solutions are contained in refer- 
ence 14. The notations in this reference are the same as the ones 
mentioned in the previous paragraph in connection with the Hartree 
solutions. In that report the porous flow velocity v w is character- 
ized by 


\= v, 


w 


p u* “ v w 



From this apul the previous expression for z± 
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The value of Zj_ is presented in reference 14 as a function of X for 
two values of P (P = 0 and p = l). Therefore, equation (57) presents 
a relation between X and A for a certain value of p. Equation (57), 
therefore, connects the parameter X which determines the shape of the 
exact velocity profile with the form parameter A determining the 
shape of the approximation. Figures 9 and 10 show a comparison between 
the exact profiles and the approximations for the neighborhood of the 
stagnation point (p = l) and for a flat plate (p = 0) . It may be 
seen that again this agreement is quite satisfactory for p = 1. For 
P = 0 the approximation is less good for profiles with a distinct 
S shape. For the largest coolant flow (characterized by the highest 
value of X or the smallest value of A) the deviation is considerable, 
although the agreement is still better than the four -term polynomial 
approx ima tion indicated by a dash-dot line. Therefore, the method pre- 
sented in this report should not be used for excessive coolant-flow rates. 


Temperature Profile 

For a solid wall (v w = 0) the approximate expression (27) for the 
temperature profile has a unique form, independent of the pressure 
gradient along the surface. For a gas with a Prandtl number of 1 and 
for constant .pressure along the surface, it is known that the tempera- 
ture profile is similar to the Blasius velocity profile. This simi- 
larity also holds for the approximation in this report. In addition, 
it is shown in reference 1 that the shape of the actual temperature 
profile is influenced only to a minor degree by a pressure gradient in 
the flow. Therefore, these approximations should be quite good for a 
solid surface. For a porous wall, with a constant pressure along its 
surface, and for a fluid with a Prandtl number of 1, the velocity and 
temperature profiles are again similar. Therefore, figure 10 shows 
also the degree of approximation for the temperature profile. 

For the flow in the neighborhood of a stagnation point and a 
fluid with a Prandtl number of 1, some temperature profiles were cal- 
culated in reference 14. Figure 11 shows the exact temperature profiles 
and the approximations presented in such a way that the thickness of 
the temperature boundary layer, which is defined in the same way as the 
displacement thickness, is the same for the accurate solution and the 
corresponding approximation . Figure 11 shows that the agreement is 
satisfactory as long as the coolant flow characterized by the value X 
or At is not too large. The approximation for the highest coolant 
flow rate (X = 3.191) is not shown since it is obvious that the exact 
curve cannot be approximated by the trigonometric functions from which 
the approximation is composed. Large coolant flows in the ranges where 
the agreement ceases to exist will probably not be used since the gain 
in cooling effect for a given increase in coolant flow is too small in 
this range. ' 
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NUMERICAL CALCULATIONS FOR TURBINE BLADES 

The method developed in the previous sections was used to calculate 
the local distributions of the heat-transfer coefficients for two tur- 
bine blade shapes for which the velocity distributions were known. The 
two blades are shown in figure 12. The line at the nose of the blade 
indicates the circumferential direction and the arrow, the direction of 
the upstream velocity. Figure 13 shows the velocity distributions 
around the blade circumferences. The velocity distributions around the 
blades were calculated using the method described in reference 17. 

The method outlined in this report can be used to calculate the 
development of the laminar -flow boundary layer in the downstream direc- 
tion up to the point of separation. In reality, at the Reynolds numbers 
encountered in gas turbine engines, the boundary layer usually becomes 
turbulent before iu reaches this point. The transition point to tur- 
bulence is determined by the pressure gradients along the blade surface, 
by the boundary -layer thickness, by the temperature distribution within 
the boundary layer, and, in addition, by the curvature of the blade 
surface and by the stream turbulence (reference 18). Not too much 
quantitative knowledge exists on the influence of all these parameters. 
However, there are indications that on turbine blades the point of 
transition is near the point at which the maximum value of the velocity 
is reached. A calculation which deals with laminar boundary layers is 
therefore useful for the region near the nose of a turbine blade. The 
points of maximum velocity are indicated in figures 14 to 16 by small 
dots . 


The first pa rt of the calculation deals with the problem in which 
a uniform coolant flow rate v w is prescribed along the blade peri- 
phery. The development of the flow and of the thermal boundary layers 
was calculated in the way described in appendix C. The results are 
presented in figure 14. Apparently both boundary layers start out with 
a finite value at the stagnation point and increase in thickness in the 
downstream direction on both the suction and the pressure surfaces. The 
boundary-layer thicknesses at any place along the blade increase with 
increasing coolant flow rate. 


In figure 15 the local heat-transfer coefficients along the blade 
surfaces which are obtained from the thermal boundary-layer thickness 
are shown. The figures show that very high local values are encountered 
at the stagnation point and that the values decrease considerably with 
increasing distance from this point. The application of transpiration 
cooling decreases the heat-transfer coefficients effectively from the 
values obtained on a solid blade surface (v w = 0) . The decrease is 

not as pronounced at the stagnation point itself a6 on the sides of the 

v w * 

blade. The value of the velocity ratio may be obtained from the 
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parameter v w presented in the diagrams by dividing it by the square 

root of the upstream Reynolds number. For turbine blades, this Reynolds 

number is near the value 10^ so that to a value v w = 1 corresponds a 

velocity ratio = 0.00316. Values for this velocity ratio which are 

Uq* 


considered m practical applications are in the neighborhood of 0.01. 

A comparison of figures 15(a) and 15(b) shows that generally the values 

of Ru at the stagnation point decrease considerably with an increas- 

/ Re 0 

ing ratio of the radius of curvature at the blade nose to the blade 


chord. The values — — are determined by the velocity gradient — , 

/Re 0 clx 

which in itself depends mainly on the dimensionless radius of curvature. 
Increasing the radius of curvature at the blade nose is therefore an 
effective means for decreasing the heat-transfer coefficient in this 
region. Downstream of the point' of transition of the boundary layer to 
turbulence the heat-transfer coefficients will increase. To give some 
indication of the magnitude of the values which may be expected in the 
turbulent region, heat-transfer coefficients were calculated at the 
location x = 0.5 under the assumption that these values are the same 
as on a flat plate at the same distance x from the leading edge, with 
a stream velocity equal to the local value found at the blade at 
x = 0.5, and with an upstream ReQ equal to 10^. Formulas derived 
by Rannie and Friedman (reference 18) were used for these calculations. 
The short horizontal lines in figure 15 indicate these values and the 
dashed part of the curve v w = 0 in figure 15(a) shows the probable 
heat-transfer values in the transition and the turbulent region. 


The surface temperature of the blade is determined by the heat- 
transfer coefficients. Equation (55) gives the connection between both 
values when heat conduction in the blade wall and heat radiation can be 
neglected. Figure 16 shows the values 1 - $ determined from equa- 
tion (55). The quantity 1 - is the ratio of the difference in blade 
surface temperature minus coolant temperature to the difference in 
effective gas temperature minus coolant temperature. The higher heat- 
transfer coefficients near the leading edge create high blade tempera- 
tures at that location. Heat conduction within the blade wall tends to 
reduce these high local temperatures. This reduction is assisted by the 
fact that the blade temperatures are especi ally low in the laminar 
region on both sides of the leading edge. Therefore, it may be that the 
blade-nose temperature is reduced to a value not higher than the blade- 
wall temperature in the turbulent part, especi ally for blades with a 
blunt nose. On the other hand, it has to be kept in mind that the heat 
conductivity values for porous materials are lower than for solid walls 
(reference 18). It has been already mentioned that temperature gradi- 
ents along the blade surface influence somewhat the local heat-transfer 
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coefficients; however, not enough information is available to take this 
effect into account. A surface temperature -which decreases in flow 
direction tends to increase the heat-transfer coefficients so that the 
decrease in surface temperature on both sides of the stagnation point 
may be somewhat less than those shown in figure 16 where heat conduction 
within the blade walls is excluded. 

A second set of calculations was made for blade I to determine that 
distribution of the local coolant flow rate v w along the blade sur- 
face which results in a constant blade tenperature. The procedure for 
such a calculation is also described in appendix C. The local coolant 
flow rates v w which are necessary m order to keep the temperature 
ratio $ (equation (55)) at the value 0.7 are shown in figure 17. The 
flow boundary layer and the thermal boundary layer determined in such 
a way are presented in figure 18. A comparison of this figure with 
figure 14(a) indicates that the development of the boundary layers is 
considerably different for both cases. The boundary-layer growth is 
smaller for the case of constant wall temperature. Figure 17 shows 
that the highest local flow rates are necessary near the stagnation 
point in order to keep the wall temperature down at that place. The 
magnitude of the coolant flow rate at the stagnation point is propor- 
tional to the square root of the velocity gradient — , which is in 

dx 

itself determined mainly by the value of the radius of curvature at the 
blade nose. The smaller this radius, the larger the velocity gradient 
and the local coolant flow rate. The flow rates decrease on both the 
pressure and the suction surfaces in the downstream direction. Down- 
stream of the transition point the boundary layer will be turbulent. 

It is known that the cooling effectiveness of a turbulent boundary layer 
is less than that of a laminar boundary layer. Correspondingly higher 
values of the coolant flow rate are necessary in the turbulent region 
in order to keep the blade temperature constant. In order to give some 
indication of the magnitude of the coolant flow rate necessary for the 
turbulent region, the values v w which result in the same value 
$ = 0.7 as in the laminar region were calculated at the location 
x = 0.5. For this calculation the formulas presented in reference 18 
were used and it was again assumed that the local $ value on the blade 
surface is the same as on a flat plate at the same distance from the 
leading edge and with a velocity equal to the local stream velocity U 
and an upstream Reynolds number ReQ equal to 10^. The flow rates v w 

obtained in this way are indicated by short horizontal lines in fig- 
ure 17. The values m the turbulent region along the blade can be 
expected not to deviate much from these values. Apparently a consider- 
ably higher coolant flow rate is necessary in the turbulent portion of 
the blade surface than m most of the laminar portion. In order to 
simplify manufacturing problems, turbine blades may be produced in such 
a way that the coolant flow rate has a constant value along the pressure 
surface and a different constant value along the suction surface. This 
then results in lower local blade temperatures in the region of laminar 
boundary layers and the heat conduction into this cooler part of the 
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blade may reduce the local temperatures at the stagnation point to 
values equal to or lower than the temperatures in the turbulent region, 
especially when the radius of curvature at the blade nose is not too 
small. For blades with a very pointed nose it may be necessary to have 
a coolant flow rate at the stagnation point approximately twice as large 
as along the sides in order to reduce the temperature at that point. 


CONCLUSIONS 

An approximate method which permits the calculation of local heat- 
transfer coefficients and surface temperatures in the laminar flow 
region around cylinders of arbitrary cross section with transpiration- 
cooled walls was developed. Velocity and temperature profiles-in the 
boundary layer were approximated by trigonometric expressions. The 
method was applied to determine local heat-transfer coefficients and 
surface temperatures in the laminar region of two transpiration-cooled 
turbine blades for a given coolant flow rate. Coolant-flow distribu- 
tions necessary for maintaining uniform blade temperatures were also 
determined. The following conclusions are made: 

1. For small coolant flow rates, the assumed trigonometric approxi- 
mations for the velocity and temperature profiles within the boundary 
layer agree well with exact solutions. For high coolant flow rates 
(which are outside the range of practical interest) the agreement 
becomes poor. It is, however, still better than the formerly used 
four-term polynomial approximation. 

2. Transpiration cooling results in a considerable reduction in 
heat-transfer coefficients on turbine blades in the laminar region. 

3. The surface temperatures of a blade with negligible heat con- 
duction are highest at the leading edge, lowest in the rest of the 
laminar region, and have intermediate values in the turbulent region. 

4. The -following variation of the coolant flow rate along the blade 
surface is necessary to keep the blade temperature constant. Highest 
local coolant flow rates are required at the stagnation point in order 
to keep the blade surface temperature down at that point. The coolant 
flow rates decrease very markedly in the downstream direction on both 
the suction and pressure surfaces in the laminar region., However, they 
increase again when the turbulent region is reached. 

5. To maintain the Bame surface temperature on both sides of the 
blade, different coolant flow rates should be used for the suction and 
pressure surfaces. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio 
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APPENDIX A 
SYMBOLS 

The following symbols are used in this report: 

c chord length of turbine blade, ft 

Cp specific heat at constant pressure, Btu/(lb)(°F) 

f. function of shape factor of velocity profile (see equations (38) 
and (40)) 

f P function of shape factor of velocity profile (see equations (39) 
and (41)) 

fj function of Pr, K, v w , and 8 (see equation (46)) 
h heat-transfer coefficient, Btu/(sec)(sq ft)(°F) 

i 

K ratio of thermal to flow boundary-layer thickness, 6^./& 

k thermal conductivity, Btu/(sec) (ft) (°F) 

L length of either side of surface of profile measured from stag- 
nation point, ft 

m exponent (called Euler number) in U = U]x m 

Eu Nusselt number, hL/k 

Pr Prandtl number, Cpp/k 

p* pressure along body, lb/sq ft 

Reg Reynolds number, UqL/v 

s distance between blades, ft 

t temperature difference between local and free -stream values, °F 
U nondimens ional stream velocity, U*/Uq 

U* stream velocity, ft/sec 

Uq upstream velocity, ft/sec 
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constant in wedge -flow velocity, U » U^x 111 

nondimens ional velocity component along surface, u */Uq 

component of velocity along surface, ft/sec 

V * / 

nondimens ional velocity component normal to surface, — a /R e n 

u 0 

velocity component normal to surface, ft/sec 
nondimens ional distance along surface, x*/L 
distance from stagnation point along surface of profile, ft 

3 

nondimens ional distance normal to surface, — /\/ReQ 


distance normal to surface, ft 

nondimens ional distance used in wedge-type flow, 


, -JL.a/2 
a/2^B V* 


nondimens ional displacement thickness used in wedge -type flow, 
&•? ITT 


-B M x 


a/^F 

shape parameter for wedge-type flow, 

shape parameter of velocity profiles (see equations (23) and 
(25)) 

shape parameter of temperature profile (see equation (28) or (48)) 


flow houndary-layer thickness, ft 

flow houndary-layer displacement thickness, ft, 


f (‘ ■ i. 


flow boundary-layer momentum or impulse thickness, ft, 

f u l 1 - n) dy 


thermal boundary-layer thickness, ft 
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e 

X 


thermal boundary-layer convection thickness, 



nondimens ional temperature ratio, t/t w 

shape parameter for exact solutions for porous walls, v w 


|i viscosity, lb/(ft)(sec) 

V kinematic viscosity, p/p, sq ft/sec 

p density, lb/cu ft 

<1> nondimens ional temperature ratio, t w /t c 


Subscripts : 

c coolant (at blade surface) 

0 upstream condition 

w wall 

Superscripts: 

* dimensional 

1 denotes differentiation with respect to x 
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APPENDIX B 


DETERMINATION OF THERMAL BOUNDARY -LAYER CONVECTION THICKNESS 


The determination of the thermal houndary-layer convection thick- 
ness 5+ which is similar to the impulse thickness of the flow 
x> ^ c 

boundary layer, results from the evaluation of equation (15) 


6 


t,c 



egdy 


(15) 


when the assumed temperature and velocity profiles are inserted in the 
integrand. A single temperature profile was assumed; it is given by 


equation (27) and is for 



6 o 1 = sin 


and ■ for > 1 



5 ip 




(27) 


J 


Distinct velocity profiles were assumed, depending on the sign of the 
shape parameter A. Moreover, these profiles were further subdivided 

y y 

into different expressions depending on g- < 1 or A > 1. These 
p- o 

assumed velocity profiles are given by equations (22) and (24) ; they 
are for A > 0, 0 <. ^ < 1 


/ - 3 Z 

u j X y a -i 5 . n y 

. e -Bin j | 


and for £ > 1 
■ o — 




u 
U = 


n = 1 - Ae 


3y 

- "5" 




( 22 ) 
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for A < 0 , 0 < | < 1 

u . « y , . . « y / . • « y\ 

o= sln sr + Asin 2 6 V 1 - Sin js) 


and for — > 1 
o 


“.I 

U 




J 


( 24 ) 


It is at once obvious that, for the evaluation of 6^. Q , 
cases must be considered. * 


four distinct 


The upper limit of the integral in the evaluation of 8^ c may be 

restricted to the value 8-fcj by virtue of the assumed temperature 
profile. If 5 -j- < 5 , the value of 8^ c can be determined from a 

single integral, with 0 and 8^ as limits, because the integrand 
vanishes for y > 8^. On the other hand, if 8 < 8^, two integrals 
are required for the evaluation of 8-fc,c because of the assumed 
velocity profiles. These integrals have limits 0 and 8, and 8 and 


Results of the evaluation of 8-^ c for the' four cases follow: 
Case I: A <0 and & t < 8 

V 


5 t,c ~ 




1 - sin % h - *b 6in % K7 I 1 - sin 


2 8 t 




T . Jty A . Jt y f . . Jt y\-l, 
[sin g 5 + A sin g 5 \^1 - sin ^ gjjdy 
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Integration and simplification lead to the following result: 

H? 2 - | (1 + A) (l - cos £ k) -f (k - | sin « k) + 

(1 + A + ^ + Mt) sin | (K + l) - sin f (K-l)] + 

"( K - 1) - imnj cos 

“l 1 - I) - 


2^2 j ^2 

( A + A \)[_ n (2K-l) + n(2K+l) " j(( 2K-1) 


cos 


(At + A^.) 


1 2 nK 

cos 


Kj" + rt(2+K) “ * ww 2" " «(2-K) 


rt( 2+K) 


COS 


cos rt l 1 + 


“t •*» ” K - “in * (*MJ - 3^ Bin , (K-l) - |] 

Case II: A < 0 and > B 

6 t,c - r ■ 8in i sin i ^ ^ - sin f §i)] * 

[sin £ | + A sin § | (l - sin | |Jj ay + 

[ X - sin f ^ ‘ \ 8 ln 5 ^ if ' 8 ln 5 ^)] (1) dy 


CO 

00 

00 

OJ 
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Integration and simplification in this case lead to 
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Integration and simplification for this case lead to 
6 t,c 


A \ 2 / -i A x 2 rtR A -3K A , 

& » (1 - A) - - (1 - A) - cos — + AK + ^ e - j + 


( -1 + A - A t + AA t ) 


(A + AA,.) 


M 


2Kit 

(361? + n 


rt(K-l) 


K 


iy Bln “ " sin i (K+1) I + 


->] 


2K 


12K & -3K 

e 


2 ) rt (361?+ j?) 


(At - A^.) ^ 2+K j + rt ( 2 _K) " *(2-K) COS 

1 / K\ 2 nKl 

n(2+K) cos jt ^ 1 + 2 J - n cos g J 


*(-§)- 


A^t 


K 


it 


2 6(9K 2 +« 2 ) 6(9K 2 + rt 2 ) 

Case IV: A > 0 and 6^. > 6 


-3K , 31? -3K 

e + o' o^ e 

■2 . JA\ 


(9K^ + 


P 6 


’t^c 




[ x - ° in l^ - ** " sin ^)] x 


slo f 2 -i- a(i - e 6 - sin | 2) 


dy + 




r>R 

& t 

- sin J JL 

A • * 

- A , sin 7 T 


, .it 

1 - sin 77 

i)l 

L 

2 5 t 

t 2 

8 t 1 

^ 2 

8 t/J 


m 
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Integration and simplification for this case result in 


A - 3 £ ^ A Z'K 1 , K , Jt \ L 
3 + ^t\2 = 2 + 2jt S3J ^ K ) + 


(-1 + A - At + ^) 


K 


Jt(K-l) 


sin 


T it' (Krl)l K _j_ f n 

|_2 K J «(K+1) Sin (_2 


(K+l) 

K 


(A + AAt) 


2Kh _ 2K 2K jt 12Kr -3K 

(36K 2 + Jt 2 ) * n C " 2K (36K 2 +jt 2 ) 6 J 


(At - AAt) 


1 , 1 1 
rt ( 2+K) Jt(2-K) ~ jt ( 2 -K) 


jt(2-K) 1 jt( 2+K) 

COS ~ 2K~ r W 008 ~ 2K - 


AAt 


1 K . Jt • Jt 2 3K 2 -3K 

2 " 2„ Bln K - 6(9K? + ‘n 2 ) (9K 2 + „ 2 ) * 


-3K 


6(9K 2 +jt?) 
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APPENDIX C 


PROCEDURE FOR CALCULATION OF FLOW AND THERMAL BOUNDARY LAYERS 

In this appendix the procedure will he explained hy which the flow 
and thermal boundary layers, the heat-transfer coefficients, and the 
temperature distribution around a body of given shape may be obtained. 
As outlined previously, the boundary conditions on the surface of the 
body may be prescribed in two different ways, either by prescribing the 
coolant velocity v w * or by prescribing the blade temperature around 
the surface. The calculation procedure becomes simpler in the first 
case and this will be considered first. In addition to the value v w *, 
the distribution of the stream velocity U* just outside the boundary 
]3ayer around the body has to be known. This velocity distribution may 
be obtained by any of the known calculation procedures for friction- 
less flow on bodies without flow separation or it has to be determined 
experimentally when flow separation occurs. From the velocity U*(x) 
the gradient U*'(x) can be found. ■ , 

The values mentioned have to be transformed into the dimensionless 
quantities U(x), U'(x), and v w with the help of equations (4). 

I 

Flow Boundary Layer 

The aim of this calculation is to determine the momentum thickness 
5^ as a function of the distance x from the stagnation point as 
measured along the surface of the body. 

As a ' first step, the value 5^ has to be found at the stagnation 
point. For this purpose determine v w ^/u ' , and read the value U'8-j^ 
from figure 4. This determines the momentum thickness 6j_. The 
intersection of U'b^ with the stagnation line in figure 3 deter- 
mines A, the shape parameter for the velocity profile. The velocity 
profile itself may be determined from equation (22) or equation (24) if 
it is desired. The required parameters may also be taken from table I. 

The momentum thickness along the profile surface is found from 
equation (37) . This first-order linear differential equation may be 
solved by any of the known procedures. The method of isoclines was 
used for the solution of the numerical examples described previously. 
The use of this method will be outlined here. Equation (37) determines 
the direction of the tangents to the different 5^ curves which 
satisfy this equation. The task is to find the one curve which con- 
tains the 5 1 value calculated previously for the stagnation point. 
Figure 19 shows the 6^ values as ordinate and the distance x from 
the stagnation point as abscissa. The directions of the tangents may 
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be obtained from equation (37) for any point in this figure character- 
ized by a pair of values x and Bj 2 and may be inserted on the ordi- 
nates through a chosen sequence of distances x 2 , x 3 . . . along 

the abscissa for several 6,- 2 values. The calculation proceeds for x 
a r id a chosen 6^, in the following way: determine U'6^ and v^S^ 

and read from figure 3 the value A and from figure 2, f^ and f 2 . 
Calculate from equation (37), d(6 i 2 )/dx. Now the direction of the 

tangent may be inserted in figure 19. The same calculation is repeated 
for other 5^2 values and the corresponding directions of the tangents 
are inserted in figure 19. The same calculation is repeated for the 
distances x 2 , x^ . . . If the distance between these values is chosen 

small enough, there can be drawn in the figure with good accuracy a 
curve which starts out with the predetermined value B^ 2 at the stagna- 
tion point and is tangent to the straight lines inserted in figure 19. 
The calculation may be shortened considerably when the curve is 
inserted step by step after the tangents have been calculated for any 
value x. Then the correct value B^ 2 for the following distance x 
may be guessed from the shape of the curve up to that point and the 
tangents need to be calculated only in the neighborhood of this value. 

The values which will be needed for the calculation may be arranged 
in a table such as the following: 


X 

U 

U 1 

v w 

»i 

A 

B 

0 

x i 

x 2 








The momentum thickness B^ in this table results from the preceding 
calculation. The shape parameter A was found in the course of this 
calculation and the value 6 may be determined from figure 2. 


Thermal Boundary Layer 

The calculation now has to be restricted to a fluid with a certain 
Prandtl number. The value which characterizes the thermal boundary 
layer and will be determined in this section is the convection boundary- 
layer thickness 8-t c* this value at the stagnation point, 

determine VyB^ from the preceding table, read from figure 7 or table II 
8^. c /& and fj. This determines &t,c* Equation (52) determines the 

local Nusselt number Nu, and for the case in 'Which the heat-balance 
equation (54) is satisfied, equation (55) determines also the tempera- 
ture ratio $. The shape factor determining the shape of the tem- 

perature profile according to equation (27) may be determined from 
figure 5 and equation (48) . 
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The convection boundary- layer thickness along the surface of the 
profile is determined by equation (46). This equation may be solved in 
the same way as the corresponding flow equation. Figure 20 indicates 
the procedure. The product U&-^ c is plotted over the distance x 

from the stagnation point. The short straight lines in the figure again 
indicate the direction of the tangents. It is known from exact solu- 
tions that the gradient of any boundary- layer thickness at the stagna- 
tion point is zero and that the boundary- layer thickness itself has a 
finite value. Thus, the curve U5-t,c in fi 6 ure 20 starts out with the 

value zero and with the inclination for x » 0. In order to 

find the direction of the tangents for any value x, for instance x-^, 
assume a value U&.^ c > calculate &.^ c and &t,c/^> rea< ^ 

value K from figure 6, determine the product Pr v^BK, read f^ 
from figure 5, and calculate from the heat flow equation (46) the 

gradient ^ (UB^ c ) . Insert the corresponding tangent into figure 20 
and repeat the calculation for other values b&-t,c* Now " ttie curve 
UB+ - which starts out with the value zero for x = 0 and with the 

u y C 

corresponding inclination may be extended to x^. The calculation is 
now repeated for the next distance Xg, the curve is extended, and 
this procedure is repeated until the curve is known for the whole 
length x. The result of this calculation gives &t,c an<1 f 3* These 

values may be inserted in a table such as the following: 


X 

6 t,c 

to 

<H 

Nu 

9 

/\j Re o 

X 1 

x 2 

X 3 



i 



The ratio Nu///ReQ which characterizes the local heat-transfer coef- 
ficient is found from equation (52). For the special case that the 
heat ba lan ce equation (54) is fulfilled, equation (55) determines the 
temperature ratio 4>. 

The calculation is more tedious when the surface temperature 
instead of the flow velocity v w is prescribed along the blade. When 
the heat ba lan ce equation (54) is satisfied, the temperature ratio $ 
may be determined as a function of x. Now the momentum equation (37) 
and the heat flow equation (46) have to be solved simultaneously by a 
trial-and-error procedure. A value v w is assumed. The calculation 
procedure described previously is carried out to determine a tempera- 
ture ratio $ corresponding to the assumed v w . The calculation has 
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to be repeated until the determined $ value matches the prescribed 
value. This calculation has to be carried out stepwise for the dis- 

> x j ) • • • 

In some cases heat may be transferred to the surface of the pro- 
file by radiation and may flow along the wall of the body by conduction. 
Then an equation which takes these processes into account has to replace 
the equation ( 54) . The determination of the coolant flow velocity v w 
which results in a prescribed blade surface temperature may be deter- 
mined in the same way as was just described. The procedure, however, 
becomes quite tedious in this case. 


tances 


* 1 ' *2 
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TABLE I - STAGNATION POINT VALUES FOR FLOW BOUNDARY LAYER 


v w^i 

A 

U'S-jS 

5 

»d 

S 

f l 

f 2 

v v 5 

-0.4999 

1.0000 

0.0001 

0.1667 

0.3333 

4.0000 

0.5000 

-2.9988 

- .4 

0.78 

-.0008 

.1614 

.3399 

4.1059 

.4336 

-2.4783 

- .3 

.6065 

.01935 

.1567 

.3452 

4.2029 

.3820 

-1.9145 

- .2 

.486 

.0345 

.1532 

.3488 

4.2768 

.3471 

-1.3055 

- .1 

.411 

.05225 

.1509 

.3510 

4.3260 

.3257 

-0.6627 

0 

.359 

.0715 

.1492 

.3526 

4.3633 

.3110 

0 

.05 

.329 

.0812 

.1484 

.3533 

4.3800 

.3050 

.3380 

.1 

.3075 

.09025 

.1475 

.3541 

4.4007 

.2966 

.6780 

.2 

.279 

.1105 

.1466 

.3550 

4.4216 

.2888 

1.3643 

.3 

.242 

.130 

.1453 

.3561 

4.4508 

.2785 

2.0647 

.4 

.225 

.151 

.1447 

.3566 

4.4644 

.2739 

2.7643 

.5 

.208 

.1718 

.1441 

.3571 

4.4781 

.2692 

2.4698 

.6 

.1835 

.19175 

.1433 

.3579 

4.4976 

.2627 

4.1870 

.7 

.172 

.213 

.1429 

.3582 

4.5066 

.2596 

4.8985 

.8 

.161 

.234 

.1425 

.3585 

4.5158 

.2566 

5.6140 

.9 

.151 

.255 

.1421 

.3588 

4.5250 

.2539 

6.3336 

1.0 

.142 

.276 

.1418 

.3591 

4.5324 

.2515 

7.0522 

1.1 

.135 

.2975 

.1416 

. .3593 

4.5374 

.2497 

7.7684 

1.2 

.128 

.3195 

.1413 

.3595 

4.5442 

.2478 

8.4926 

1.3 

.123 

.341 

.1411 

.3597 

4.5492 

.2464 

9.2133 

1.4 

.114 

.362 

.1408 

. 3599 

4.5561 

.2449 

9.9432 

1.5 

.113 

.384 

.1408 

. 3600 

4.5568 

.2439 

10. 6534 

1.6 

.107 

.405 

.1406 

. 3602 

4.5618 

.2424 

11.3798 

1.7 

.101 

.426 

.1403 

.3604 

4.5688 

.2406 

12.1169 
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Pr K v w 8 



Figure 5. - Chart for determination of f 3 used in thermal 
boundary-layer differential equation. 
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Figure 6. - Chart for determination of K for thermal boundary- layer calculations, (A 17 « 
by 20- in. working chart of this figure is enclosed.) 
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Figure 7, - 


Chart for determination of 6 t>c< & and f 3 at stagnation point for thermal boundary-layer calculations. (A 17- 

by 9- in. working chart of this figure is enclosed.) 
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Figure 8. - Comparison of approximate velocity profiles for solid surface with Hartree solutions (reference 16). v w , 0. 
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Figure 10. - Comparison of approximate velocity profiles for flat plate with Schlichting solutions 
. ( reference 14). fl, 0. 
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Figure 11. - Comparison of approximate temperature profiles for flow at stagnation point with 

exact solutions (reference 14) 
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(a) Blade I* Ratio of chord length to distance between consecutive 

blades c/s, 1.408. 



(b) Blade II. Ratio of chord length to distance between consecutive 

blades c/s, 1.136 


Figure 12. - Turbine-blade profiles used for calculated examples. 
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' (b) Blade II. 

Figure 16. - Surface temperature distribution for blade. Pr, 0.7 
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Figure 17, - Variation of coolant flow required to maintain constant 
blade temperature for blade I. 0,7, Pr, 0.7. 



Figure 18. - Growth of flow impulse and thermal convection 
boundary-layer thicknesses for constant blade temperature 
for blade I. 0, 0.7; Pr, 0.7. 
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